arXiv:hep-ph/0503013vl 1 Mar 2005 


Speed-up of neutrino transformations in a supernova environment. 
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When the neutral current neutrino-neutrino interaction is treated completely, rather than as an 
interaction among angle-averaged distributions, or as a set of flavor-diagonal effective potentials, 
the result can be flavor mixing at a speed orders of magnitude faster than that one would anticipate 
from the measured neutrino oscillation parameters. It is possible that the energy spectra of the 
three active species of neutrinos emerging from a supernova are nearly identical. 
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INTRODUCTION 

Predictions regarding the neutrino flavor physics in 
the region of the supernova’s neutrino-sphere have re¬ 
mained quite stable as treatments of other issues have 
evolved over time. These predictions, which involve sig¬ 
nificantly different energy spectra for the u e , v e , 

P/r,r) are important to the theory of the region above the 
neutrino-sphere. They matter for the explosion dynam¬ 
ics, for heavy element nucleosynthesis, and for analyzing 
the neutrino burst from the galactic supernova scheduled 
for some time in the twenty-second century. 

As the three known species of neutrinos move through 
layers of the supernova that are immediately under the ef¬ 
fective neutrino-sphere, their physics can be much more 
complex than is assumed in current simulations of the 
whole supernova process. Before describing some of these 
complexities in detail, we look at parameters that de¬ 
scribe the region: taking a density of 10 11 gc -3 and tem¬ 
perature of 7 MeV. There is a characteristic energy that 
we define from the Fermi constant and the neutrino num¬ 
ber density, E s = 2^/2 Gfti 1 ,, where n v is the neutrino 
number density, which when translated into distance, 
under the above conditions, gives Ej 1 = ,04cm 1 . Of 
course, this is many orders of magnitude less than the dis¬ 
tance set by the density and cross-section (proportional 
to Gp 2 ) , a distance scale that is tens of km., the scale 
of the neutrino-sphere radius itself. 

However, terms of order E s do enter in forward (or “in¬ 
dex of refraction”) effects. In a medium containing elec¬ 
trons, the v-e interaction energies create different poten¬ 
tials for i/ e ’s, and for P e ’s, than for and v T 's and their 
antiparticles. The induced terms of order E s enter inter¬ 
esting physics in conjunction with the (5m 2 ) neutrino os¬ 
cillation parameters, but these oscillation effects do not 
enter significantly out to a distance of many times the 
neutrino-sphere radius, if we use current values for the 
neutrino mass parameters and operate in the commonly 
accepted theoretical framework. The electron densities 


1 We use units in which h = 1, c = 1 


have no effect on the oscillation between fi u and /i T , with 
oscillation length of the order of a kilometer for a 20 MeV 
neutrino, and uninteresting in any case by virtue of the 
near identities of their spectra snd angular distributions. 
The v e oscillation length in vacuum is longer by a factor 
of 40 or so; furthermore, the large electron density en¬ 
sures that oscillation amplitudes will be miniscule. Thus 
there is negligible “forward” neutrino physics of any con¬ 
sequence in the conventional picture of this region, bar¬ 
ring a sterile neutrino that is attached to a much bigger 
mixing parameter. 

All of this can change when we include the full neutral 
current interaction among the neutrinos, however. The 
effective forward interaction Hamiltonian for the system 
now contains terms that can exchange neutrino flavors. 
At first sight this seems fairly irrelevant in the region 
below the neutrino-sphere. The initial system coming 
from the supernova core can be taken as an incoherent 
sum of flavor eigenstates. Then the mixing of the v e re¬ 
mains small, according to the usual angle-averaged equa¬ 
tions, which assume isotropic momentum distributions. 
A large, nearly flavor diagonal, surrounding neutrino den¬ 
sity can serve to synchronize [ij-js] these small oscilla¬ 
tions, according to these equations, but not to enhance 
their amplitude when the amplitudes otherwise would be 
small. 

In a region well above the neutrino-sphere the oscilla¬ 
tion parameters can produce strong mixing and possible 
MSW transformations. Considerable attention has been 
given to the effects of the neutralcurrent coupling terms, 
includin g th e forward flavor-exchange parts, in this re¬ 
gion |9(-[l2j. Here neutrino anisotropy was taken into 
account by using the isotropic form for the equations, 
but with an effective coupling constant that incorporates 
a single (“flux-averaged”) angle of incidence for v-v en¬ 
counters. However we now know that other non-linear 
phenomena can take over when anisotropy is incorpo¬ 
rated into the evolution equations in a more correct way 
H3 These phenomena, to some degree, can see the time 
scale t s = Ef 1 and overturn our intuitions about the 
problem. 

In the present paper we look at the neutrino physics 
only over short distances, over which the surrounding 
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density and temperature can be taken to be constant. 
Even this abstracted problem is extraordinarily complex. 
The simulations in Jl3j|, backed up by analytic consid¬ 
erations, were based on two-flavor models with simple 
anisotropies, essentially a few clusters of definite direc¬ 
tion in momentum space. They show the possibility of 
one kind of instability, in the form of an exponentially 
increasing growth in flavor-mixing amplitudes seeded by 
the tiniest (coherent) initial mixing (which could come 
from earlier neutrino oscillation effects). In these exam¬ 
ples the rise-time of the exponential is the shortest time 
scale in the problem, of order Ej 1 . This behavior is 
in contrast to the oscillatory behavior in the isotropic 
approximation. That said, in the domains considered 
in ref. 0 the conditions that can lead to this runaway 
evolution are delicately dependent on the initial flavor- 
momentum correlations, and are subject to change when 
new elements are added to the model. 

Approaching an authoritative answer in a more real¬ 
istic supernova environment requires the inclusion of all 
three active flavors and their antiparticles. 2 Let us think 
of the neutrino distributions in the transitional region 
just below the average depth of the last v p and v T scat¬ 
terings, a region that is somewhat farther below the aver¬ 
age depth of the last v e scattering. In this region the v e 
momentum distribution is fairly isotropic, whereas the 
Pfj,, v T distributions are already quite biased to¬ 
ward upward momenta. Both in angular distribution 
and in spectrum, the P e ’s lie between the v e 's and the 
z/ MiT ’s. In our first idealization we picture this mixture 
as a core group of all six species in equal number with a 
momentum distribution peaked somewhat upwards, su¬ 
perimposed on a group in which we have ’s and their 
anti-particles in a distribution peaked strongly upward, 
and i/ e , v e 's in a distribution peaked downward. We call 
the second group the valence group. The core group is 
irrelevant to the flavor evolution of the system, as may 
be obvious but in any case is supported by our later for¬ 
malism. The problem that we now solve is one in which 
we replace the valence system by narrow angle beams of 
and v T neutrinos and antineutrinos directed outward, 
and beams of i/ e ’s and u e 's directed inwards. In a second 
version we shall assign the P e ’s to the outgoing group, 
corresponding to a less outwardly biased core group. 

Over the time-scales that we consider below, the in¬ 
teractions that come into play, in addition to the neutral 
current couplings, are the neutrino mass-mixing terms 
and the difference between the forward e-u e and e-Vp tT 
scattering amplitudes. These interactions would not even 
matter over our time scales, but for the instabilities in 
the non-linear equations that govern the evolution of the 


2 Inclusion of a sterile species would complicate our considerations 
enormously. 


system. The first initial state that we pick consists of 
down-moving i/ e ’s and P e ’s, each with a number density 
of n v , with the other four species of particles moving up¬ 
ward each with the same number density. The outcome, 
in this case and in two other cases with somewhat dif¬ 
ferent initial conditions, is a scenario of extremely rapid 
mixing of flavors, leading to a rapid destruction of some 
of the initial correlations of flavor with angle, and in turn 
with energy spectrum. 

EQUATIONS OF EVOLUTION 

The equations are an extension of the equations of Sigl 
and Raffelt 0- The main new element is the explicit 
inclusion of all three flavors of active neutrinos and an¬ 
tineutrinos, as well as perhaps a more incisive inclusion of 
angular factors, so that we can better treat anisotropy ef¬ 
fects. The formalism is much simpler than that presented 
in ref. 0 because of our limitation to purely forward ef¬ 
fects, justified by the short time scales to be considered. 

All of the dynamics of the system can be described 
in terms of the annihilation and creation operators for 
the three neutrino flavors and the momentum states of 
the initial configuration; we define the annihilators of 
the mode p for the respective three flavors u e , v T , as 
a p ,b p ,c p with the annihilators for the antiparticles de¬ 
noted ap,bp,c p . The neutral-current effective interaction 
Hamiltonian, to be exhibited below, is independent of the 
neutrino energies and depends only on relative directions 
of motion. Therefore, for the purpose of discussing the 
evolution of the density matrix due to this interaction 
alone, we can bundle all energies together and define col¬ 
lective operators for each angle. We divide the solid angle 
of the momentum space into elements dfl centered at an¬ 
gle Q and define collective operators for each of these 
elements using a shorthand symbol Sq which operates 
according to the rule, 

5n[a t 6] = (dQ) _1 ^ aj,6 p . (1) 

p C dQ, 

We choose, 

pi(Q) = S'o[a t 6] , P2(H) = Sn^a] , 
p 3 (fi) = S'n[(2c t c - a)a - tfb)/3] , 
p 4 (fl) = , p 5 (Q) = 5fi[c f 6] , 

Pe{^) = S’n[(2a t a - tfb - c^c)/3], 

MH) = <S’n[c t a] , ps(^) = , 
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p 9 (f2) = S'n[(26 t 6 - at a - c*c)/3] , 

Po(^) = S'nta^a + b^b + c^c]. (2) 

Given the redundancy P 3 + pe + P 9 = 0, these embody 
the usual SU3 relation, 3* x 3 = 1 + 8. We construct 
bilinears of the antiparticle operators a, b , c in the same 
fashion as above except that we make the exchanges in 
indices, relative to the neutrino case, of 1 <-> 2, 4 <-> 5, 
and 7 <-> 8. For example, we take, 

pi(Cl) = Sn[P a] , p 2 (fi) = Sn[at b] , etc. (3) 
The commutation rules are now given by, 

9 

\pi(Q),pj({}')} = 5(SI - SI') ^2 fij,kpk(fy , 

k =l 

9 

[pi(fi),Pj(n')] = - n ; ) y; fij,kpk(s}) , 

/e=l 

[p < (n),p i (n , )] = o. (4) 

The fij.k coefficients characterize the SU3 algebra of 
the operator arguments of the expressions m as they 
describe a single momentum state. They are listed in the 
appendix. 

In terms of these operators the forward neutral current 
interaction Hamiltonian is, 3 

H ior= y^voi.) f dn dV'il - cos6> n ,o') 

{§[po(f2) - pom x M^') - Po(O')] + 

E hlPim - Pi(n)]b4(n') - pt(n')]}. (5) 

The result © is obtained from the full neutral current 
four-Fermion interaction by retaining only those com¬ 
binations of operators that give a forward amplitude, 
p, q —> p, q , and restricting these momenta to the set 
of states that are initially occupied. The terms in © for 
v-v (and v-v) scattering come from two graphs. In one 
of these a zero-momentum Z is exchanged, giving rise 
to an interaction among the po and po operators only, 
since for this term in the scattering p + q —> p + q, the 
operator that makes the final p and the operator that 
annihilates the initial p attach to the same end of the 
Z line. This singlet-singlet coupling (in the SU3 sense) 
gives 3/4 of the first term in © and contributes noth¬ 
ing to the time evolution of the density matrix. In the 
other (“crossed”) term for v-v scattering, a momentum 
p — q is transferred through the Z. In this term we use a 


3 The kinetic term in the Hamiltonian takes the same value for 
every state in our set, and therefore can be ignored in all that 
follows. 


Fierz identity, giving a multiplying minus sign, one anti¬ 
commutation, giving a second minus sign, and finally we 
use the appropriate crossing matrix in the internal SU3 
space to produce the rest of the terms in ©. For the 
particle-anti-particle terms, the crossing is slightly dif¬ 
ferent, since the “crossed” graph is a virtual annihilation 
into an intermediate Z, in which the Z carries momen¬ 
tum p + q. Our 1 <-> 2, etc., exchange in subscripts in 
defining the anti-particle densities reflects this change. 
Of course, since the neutrinos will be in the MeV energy 
range, we can neglect the neutrino momenta in the Z 
propagator. 

To the above Hamiltonian we add a neutrino mass 
term, taking the v and the v's all to have nearly enough 
a common energy, E, so that the oscillation term is the 
same for each. (We weaken this assumption later.) We 
choose, 

i^mass = W / dfl[pi(fl) + + Pl(^) + P 2 (^)] + 

A 2 / dfi[p4(n) + P5(fl) + P4:(Q) + P5(fl)] , (6) 

where Ai is given by (2 E)~ x times the 5m 2 parameter 
for v e T mixing and A 2 is given by (2 E)~ x times a 5m 2 
parameter for v^^ T mixing. Finally we include the (rela¬ 
tive) energies of interaction that result from a net n e — n e 
electron-positron density difference, 

H e = \f2GF{n e - n e ) J d^l[pe{fl) - pe(fl)}. (7) 

From the above Hamiltonian and commutation rules 
we can write Heisenberg equations for our set of opera¬ 
tors Pi{ Cl), Pi(fl) . The right hand sides will be quadratic 
functions of the densities themselves, with one angular in¬ 
tegration. Thus if we are allowed to replace the expecta¬ 
tion values of products of density operators, on the right 
hand sides, by products of expectation values, we obtain 
a closed set of non-linear, integro-differential equations 
for the expectations. More explicitly, we take an initial 
condition in which the expectation of any product of our 
density operators is the product of the expectations of 
the individual operators, e. g., 

{p 3 {fl),t)p 6 {n,t)) = (p 3 {n,t)){p6{n,t )), (8) 

at t = 0. The further assumption is that this factor¬ 
ization holds for all subsequent times. In the present 
work we assume that this is justified. 4 To illustrate the 


4 Truncation schemes such as this should be regarded with some 
suspicion although we note that the whole literature pertaining 
to the effects of v — v interactions is ba sed on the factorization 
assumption. In previous work m, m, in cases of simpler al¬ 
gebras, we have found that the defect of the expect at ion-value 
assumption can be expressed as additional terms proportional 
to ATy 1 in non-linear equations for the products of density op- 






complexity of our system, we write the first of the ten in¬ 
dependent equations of evolution, where now all density 
variables are understood as expectation values, 
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Vi(fi) = (1 - cos0n,n'){p?(^)[^(^) - 

Pim + pi(tt)Mfr) - P6 (n') - p 9 (n /) + P6 (n')] - 
p 4 (nm(n') - pW)} + \pein) - P9 (n)][ Pl (fi') - 

)oi(n')]} + Ai [P6 (^) — Pg(f 2 )] + Ag/s^fl) — 

\f2GF{n e - n e )pi(Q). 

(9) 

The equations are unmanageable except in simplified 
situations. As mentioned above, there have been a num¬ 
ber of treatments of the two-flavor case in which the an¬ 
gular distributions are isotropic. In work reported below 
we have found non-linear phenomena for the three-flavor 
case that are qualitatively different from the “synchro¬ 
nized” oscillations that can arise in the two-flavor sys¬ 
tem, even in the isotropic case. We return to this point, 
but we include anisotropy in all of the detailed solutions 
that we present later. Now we implement the two-beam 
idealization, as discussed in the introduction. Setting 
1 + cos0n,n' = 2, we define collective density coordinates 
for the two beams as, 


TT 2V2N2G Ff 4 I 

H i ° r = ( Vol ) —IgR) - x o\[yo - yo\ 


+ Y^[xi -XiM -2/1]} • 


( 12 ) 


The i/-mass and electron interaction terms are given by, 

-ffmass + He = N v 2\f2GFn v [\i{Xi + X2 + 3/1 + 2/2 + 

xi + x 2 + yi+ §2) + + £5 + 2/4 + 2/5 + £4 + 

x 5 + 2/4 + 2 / 5 ) + {n 2 n™ e) (^6 + 2/6 - Xq - y 6 )} , ( 13 ) 

where dimensionless oscillation parameters have been in¬ 
troduced, A' = Xi /(2v / 2G Fn „). Writing the Heisen¬ 
berg equations that come from the total Hamiltonian 
Hfor + -ffmass + H ei and the commutation rules II fill and 
measuring time in units (2\/2G fti v )~ 1 , we obtain a set 
of 24 coupled nonlinear equations for the 24 density ma¬ 
trix components. We show the first two here; more are 
provided in the appendix, 


ix 1 = £1(2/9 - 2/6 - 2/9 + Ve ) ~ £4(2/7 - 2/7) + 
£7(2/4 - 2/4) + (£6 - £9)(2/1 - 3/1) + K ( x 6 - £9) + 


^2 X 7 


n e —n e 
2 n v 


£1 , 


Xi = {N v ) 1 / dfl pi(fl ), 

J down 

2 h = (N„)~ 1 f dClpi(Cl ), ( 10 ) 

J up 

where N v = n,,(Vol.) is the initial number of each flavor 
of iz, in our “valence” group, and n v is the corresponding 
number density. We define x^, yi in parallel fashion. The 
commutation rules are now of the form, 

[X {, Xj ] (A^jy) ^ ( fi,j,k x k 7 

k 

\Xi 1 Xj ] (-^ 1 /) ^ ^ fi,j,k%k •> (11) 

k 

The yi , yi variables have the same commutation rules 
among themselves, and \x^yj\ = 0. 

The operative neutral-current Hamiltonian for these 
variables is, 


erators. However we also found that when the non-linear back¬ 
ground equations had a certain instability, these terms could lead 
to rates of otherwise forbidden processes that are suppressed only 
by a factor log(AT), and may need to be retained even when N 
is huge. 


a 3 = £4(2/4 - 2/4) - £4(2/4 - y %) + 

£7(2/7 - 2 / 7 ) - £7(2/7 - 2/7) + A' 2 (£| - Xi ) . (14) 

We generate the non-linear parts (i. e., the parts com¬ 
ing from neutral current interactions) of the 22 remain¬ 
ing equations (four of them redundant) from the following 
successive steps: 1) simultaneous permutations of indices 
123 —> 456, and 456 —> 789; 2) the same step repeated; 
3) the interchange x <-> y (bringing us up to a total of 
12 equations) ; 4) adding a bar to each unbarred symbol 
and removing the bar from each barred symbol, com¬ 
bined with an all-over change in sign of the right hand 
sides. The remaining linear terms on the right hand side, 
containing the I'-mass terms and the electron interaction 
term are better computed by hand from m using CD- 

These equations are the gen eralization of those derived 

in ref. 13, and applied in 1-13 , for a situation which 

has only two flavors and which is spherically symmetri¬ 
cal. The equations of ref Q, specialized to two flavors 
and to a single energy, can be easily recaptured from ©■ 
We can improve a little on the set m ■ as extended by 
the permutations listed above, by remembering that the 
up-moving particles have, on the average, higher energies 
than the down-moving ones. We incorporate this through 
reducing the oscillation terms for the y variables by a fac¬ 
tor £ = (E e )/(E IJ-jT ), as calculated using the initial values 
of these energy averages, leading to equations of the form, 
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iyi = Vi(xg - x 6 - Xg + x 6 ) - 2/1(zy - x$) + 
y 7 (x% - x%) + (ye - yg)(xi - ah) + ^X[(ye - yg) + 

^ A 2^7 2n v y 1 ’ 

iy 3 = y%(xi - xa) - y 4 (a ;2 - x|) + 2/7(37 - *7) 

~t,y 7 (x 7 - x 7 ) + £\' 2 (yl - Va) , etc. (15) 

We take £ = .6 in what follows. Turning to the solu¬ 
tions, we take a neutrino density for each species, n v = 
18(MeV) 3 , about 1/2 the thermal density for a temper¬ 
ature of 7 MeV. Thus our picture is that roughly half of 
the neutrinos fall into the “valence” group. Then the unit 
of energy in our scaled equations is E s = 2 >/ 2 ti v Gf ~ 
5 x 10“ lo MeV (corresponding to a distance parameter 
of ss .04 cm.) We take an average v energy of 17 MeV. 
The v e oscillation parameter, A) = dm 2 _ r /2i%, in units 
of E a , is now A'j ~ 6 x 10“ 9 , for <5m 2 _ r = 10” 4 eV, a 
number quite unfriendly to computations since there are 
oscillations in the key flavor-changing densities at periods 
of order Eg 1 . We take X' 2 = 50Ai- For a mass density of 
8 x 10 lo gc“ 3 and with Y e = .4, we have n e /n v w 8. 

In fig. 1 we show a solution to the equations with ar¬ 
tificially large values of the oscillation parameters, A) = 
10“ 5 , and A 2 scaled to A) as described above. We take 
initial conditions that correspond to six groups of equal 
density, with u e , v e moving downward and u^ tT , mov¬ 
ing upward. This initial condition translates into initial 
values, X 3 = X 3 = xg = xg = —l/3;xe = xg = 2/3; and 
2/3 = 2/3 = 2/9 = 2/9 = 1/3; ye = y& = -2/3, with all other 
initial values equal to zero. 

From the plots of fig. 1 we see that strong mixing sets in 
shortly after a time t = (A^) -1 / 2 , with the down-moving 
states at some moments nearly totally occupied with ;/ /J T . 
Not exactly in synch with the down-moving mixing, the 
up-going states become up to 50% occupied with u e ’s, 
although the percentage fluctuates wildly with time. We 
observe that the total number of u e 's and g/s is not con¬ 
stant in the above, although their difference is constant. 
Owing to the perfectly symmetrical initial condition, the 
plots for anti-neutrinos are identical. These results are a 
prototype for all that follows; for smaller values of A), the 
time required for big mixing of the down-moving states 
follows the t = (A))~' 1 rule fairly closely as the oscilla¬ 
tion parameters are reduced. 

In fig 2. we show the case A' x = 10“' (still 16 times our 
target value). The flipping in the up-moving states oc¬ 
curs after roughly 5000/(27 t) periods of the rapid oscilla¬ 
tion rate 5 2\f2n v G f. Comparing fig. 1 with fig.2, where 


5 Luckily for our plots, this very rapid oscillation surfaces only in 
the expectations of the flavor-off-diagonal density operators. 



FIG. 1: Electron neutrino density in the up-moving beam 
(light curve) and in the down-moving beam (heavy curve) 
as a function of time, as expressed as a fraction P e of the 
original density n„. The unit of time is the fast scale, as 
defined in text, t s = (2\/2In this calculation the 
oscillation parameters have been taken to be 1.6 x 10 3 times as 
big as present oscillation parameters demand, with A) = 10 —5 , 
A 2 = 50A), in order to better show features of the curve and 
to exhibit the scaling, discussed in text, which roughly relates 
these curves to those of fig. 2. 

there is a factor of 100 difference in the oscillation pa¬ 
rameters, and a factor of 10 difference in the time spans, 
we see good confirmation of the t = (A^) -1 / 2 condition 
for mixing. Therefore it appears that we can extrapolate 
reliably to the physical value of A'. Doing so, we obtain a 
time of approximately 2.4 x 10 4 in our fast units, giving 
a mixing distance of around 1000 cm. 


Pc 



FIG. 2: The same as fig. 1, except that each oscillation pa¬ 
rameter is reduced by a factor of 100, so that A) = 10~ 7 , etc. 
The unit of time is again the short time t s , which is the oscilla¬ 
tion period of the flavor-non-diagonal operators. Comparing 
with fig. 1 we see the (A , ) 1,/2 scaling discussed in text. 

In the calculation leading to fig. 3 we take a five per¬ 
cent v e — v e surplus in the initial state, again for the case 
A'j = 10“'. We see from fig. 3 that in this case one-half 
of the up-going states become occupied with u e , D e , as in 
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the previous cases, but more steadily, whereas now the 
occupancy of the down-moving states is changed hardly 
at all. There is only a tiny difference in the behavior of 
the neutrino and antineutrino plots in this case, despite 
the fact that major features are driven by the small initial 
asymmetry. 


Pc 



FIG. 3: The same conditions as fig. 2, except for the addition 
a 5% excess of u e 's over the other species in the initial up- 
moving beam. 

In a third situation we change the rules so that the 
Ve's are placed in the upward beam. Since in typical 
supernova calculations the u e are intermediate between 
the F /Ji7 -’s and the F e ’s, both in energy spectrum and in 
angular distribution, this is a calculation complementary 
to the ones presented above. As mentioned earlier, it 
corresponds to a different choice for the “core” states. 
Again using values of A^ and \' 2 that are 16 times too 
large, because of computational constraints, we show in 
fig. 4 results for the electron neutrino distributions, in 
this case appearing to be a clean flip. In contrast to 
the earlier cases, shown in figs. 1-3, the anti-neutrino 
distributions are frozen at their initial values in this case. 


Pc 



FIG. 4: The same parameters as for fig. 2, except that in 
the initial conditions the P e ’s have been put in the up-moving 
group. 


DISCUSSION 

The purposes of this paper were: 

1) To give a simple formulation of the complete 3-flavor 
equations that govern neutrino evolution, over time peri¬ 
ods that are short compared to scattering times, when the 
density of neutrinos and antineutrinos is large. This is 
accomplished in the construction of the forward neutral- 
current Hamiltonian m and the commutation rules 0 , 
used in conjunction with the usual oscillation terms. 
Of course, in a complete approach with radically non- 
monoenergetic distributions, the densities have still an¬ 
other index, E, and we would write equations for den¬ 
sities p( E \C£). In this case the forward neutral current 
Hamiltonian has a double sum over energies but no ex¬ 
plicit further energy dependence. The energy enters ex¬ 
plicitly in the oscillation terms. 

2) To find a way of treating these equations in the 
quasi-monoenergetic case, a case that at least can give 
some insight into possible behaviors. It appears to us 
that the approximation of using a few rays in the space 
of angles of the momentum distribution, each with its 
own complete set of 20 densities is both conceptually 
clear and practical to implement, the only caveat be¬ 
ing that going beyond the two-ray approximation of the 
present paper will require large computing resources, ow¬ 
ing to the extreme values of dimensionless parameters. A 
possible alternative would be to make a partial wave ex¬ 
pansion of the density operators. Of course, because of 
the non-linearity, and the angular factor 1 — cos 00 , 0 ' in 
the kernel of ©, multipoles of all orders are coupled to¬ 
gether. Thus, even though all reactions are forward, so 
that the manifold of occupied momentum states remains 
unchanged, and even though we might choose an initial 
configuration with, say, only l = 0,1,2 states present 
in the distributions, the flavor configuration could, over 
time, become arbitrarily complex in angle, developing 
modes of higher and higher multipolarity. In any case, 
because we cannot verify the stability of the truncation 
procedure we do not pursue the partial wave approach 
here, preferring instead the discrete rays approach. In 
the rays approach, we represent continuous distributions 
in the initial state less well, but we really solve the re¬ 
sulting dynamical evolution problem. On the basis of 
the previously investigated two-flavor models we expect 
that when we begin with more than two rays there will 
be instabilities that grow on a time scale of several times 
(GfUj,) - 1 rather than {Gfti v Ai) - 1 / 2 , under the condi¬ 
tions of this present paper, perhaps rendering all of our 
present results academic 6 . Unfortunately, we are not 
equipped to deal computationally with the generaliza- 


This is also true with a truncated partial wave expansion that 
includes l = 0,1 and 2 states. 
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tion of the SU3 case, with antineutrinos included, to a 
case with more than the two rays of the present paper. 

Turning to the specific results of the present paper: In 
three different idealizations, each capturing some of the 
features that could be important in the supernova neutri- 
nosphere problem, we find an array of possible behaviors. 
In each case the time-scale is of the order (.E^Ai) -1 / 2 , 
corresponding to a reaction distance of order 100 cm. In 
the first case we took an initial condition in which both 
v e and an equal number of z7 e from the valence group 
were assigned to the down-going beam, while u T , Tv, i/^, 
v fl were assigned to the up-going beam in equal num¬ 
bers. The result was a nearly complete flavor over-turn, 
in which at some moments in time one-half of the up- 
moving beam is composed of v e 's and the down-moving 
beam is entirely composed of z^ t . T ’s. 

In the next example, we added a small extra number 
of v e 's in the initial down-moving group, uncompensated 
by anti-particles. In this case the system does a very 
amusing thing. The down-moving states remain occupied 
by v e , v e ; but the up-moving states nevertheless make 
excursions into a configuration in which half of the total 
original number of iz^’s, v T 's have transformed into z/ e ’s, 
and similarly for the antiparticles. 

In the final example we put the initial Tv's into the 
up-moving group. In this case there is perfect flavor neu¬ 
trality in the anti-neutrino sector (i.e. a nonvanishing ex¬ 
pectation only for the operator yo), and no antineutrino 
action. However, there is big mixing of the neutrinos in 
the (moderately) fast time scale, (i^Ai) -1 / 2 , as we saw 
in fig.4. 

One thing in common in all of these simulations is that 
they require both parameters Ai and A 2 to be non-zero, 
even though when Ai << A 2 << GftIv the effective 
time constant, (i? s Ai) -1 / 2 , does not depend greatly on 
A 2 . Thus all of the phenomena manifesting themselves 
here require having three flavors and two mixings. 

We can try to separate the effects of the two innova¬ 
tions in this paper: 1) the two-group treatment of the 
angular effects; 2) the inclusion of all three flavors and 
their antiparticles. We perform all of the calculations 
described above but this time beginning from the com¬ 
pletely isotropic equations. There are nearly twice as 
many terms on the right hand sides of the equations, 
as compared to 03 , owing to the intra-group interac¬ 
tions 7 In the analogs of the calculations with completely 
symmetric conditions that led to figs 1,2 , we find much 
diminished mixing effects at approximately ten times the 
time as in the case of our basic calculation. For the other 


7 A group now consists of a set of neutrinos with momenta in all 
directions. In the neutral current interaction cos 6 is replaced by 
its average value of zero. We still use two groups of states, to 
take into account the different energy spectra. 


two examples, which led to the plots of fig. 3 and fig. 4., 
the isotropic equations gave no appreciable mixing. 

It would be rash to extrapolate from any of the above 
results to draw firm conclusions for supernovas. Since 
our exploratory work 0 in the two-flavor case indicates 
the existence of instabilities that can lead to macroscopic 
changes on a shorter time scale than (E s Ai)^ 1 / 2 , once we 
introduce more complex angular dependence, someone 
with massive computing power should explore this case. 
The results could render much of the above obsolete, but 
not the basic equations or the method of attack. A reader 
could ask, “Why spend so much effort on a problem, if 
the author believes that it misses the main points of the 
ultimate solution?” The answer is that this appears to be 
quite a hard problem, and that systematic work around 
the edges may be likelier to bring us an understanding of 
what needs to be done than is plunging ahead, discarding 
this and that and approximating the other, in order to 
claim a conclusion for a physical system. 

The main moral of our present results, as well as those 
presented in ref. |j]|, is that there are big nonlinear ef¬ 
fects latent in the neutral-current couplings, effects which 
can lead to what we have called speed-ups. We speculate 
that there is a connection between the speed-ups of the 
present paper, which we would say are due to marginal 
instabilities, and those more rapidly growing instabilities 
described in ref. 0 . We think that the examples con¬ 
sidered above sit exactly on the boundary of an unstable 
region in a parameter space. Indeed, the route to finding 
the growing modes in ref. 0 was to take seriously a 
demonstration that the two-beam case (as defined in the 
present paper) is on such a boundary. If the author were 
to guess at the ultimate outcome, he would guess that 
the end result of the non-linear effects will be total flavor 
mixing, and on a distance scale shorter than the 1000cm. 
that the calculations of this paper indicate. This would 
be perhaps a rather dull outcome, but it would be compu¬ 
tationally convenient as far as incorporating the answers 
into big simulations of the complete supernova process 
are concerned. 

We emphasize that there is no MSW transformation 
in any of our results. We deal in time spans that are so 
short that the density of the surrounding medium can be 
taken as constant, in the first place. In the second place, 
the effective time for the actual turn-over, after the wait¬ 
ing time (-E s Ai) -1 / 2 , during which the system prepares 
itself, is many orders of magnitude faster than an MSW 
transition time can be for any time-changing (i.e. space¬ 
changing) medium, for the oscillation parameters that we 
choose. This distinguishes our considerations from those 
of refs. @-0]- 







APPENDIX 


( 17 ) 


£4(2/1 - yl ) + ( x 3 - xe )( y 7 - 2/7) 

AW-A^ + ^r- 

The coefficients fi,j,k which in our representation are 
antisymmetric only in the first two indices, have the fol¬ 
lowing non-vanishing components, _ 


/1,2,6 — /1,4,8 — /1,9,1 — /2,6,2 — /2,8,4 — f 3 ,5,5 — 
f 3 ,7,7 = / 4,5,9 = fi,7,2 = /5,9,5 = /6,8,8 = /7,8,3 = 1 ; 

/l,2,9 = /1,6,1 = /l,7,5 = /2,5,7 = /2,9,2 = /3,4,4 = 
/3,8,8 = /4,5,3 = /4,9,4 = /5,8,1 = /6,7,7 = /7,8,6 = —1 , 

(16) 

as well as the components obtained by interchange of the 
first two indices. Although this set is slightly less eco¬ 
nomical than the canonical set of Gell-Mann (which are 
based on a hermitean set of pi ), the present representa¬ 
tion serves us well for the purpose of writing down the 
Heisenberg equations. In this representation it suffices 
to write only the equations for the time derivatives of 
pi and p 3 , with the rest found by permuting indices, as 
we discussed earlier. The non-hermitean operators that 
we use also made it easy for us to express the quantities 
(P 2 ), (ps), (ps) in terms of the respective complex conju¬ 
gates of (pi), {pi)^ (P 7 ), reducing the number of integro- 
differential equations from sixteen to ten, or the number 
of independent ordinary DE’s from thirty two to twenty, 
in the case of our two-beam approach. 

We show the first five of these equations, beginning 
with the two presented previously. The remaining 15, 
now including the mass terms, are obtained from the final 
two substitution rules given in text: a.) the interchange 
x <-> y, b.) adding a bar to each unbarred symbol and 
removing the bar from each barred symbol, combined 
with an all-over change in sign of the right hand sides, 

ix 1 = £1(2/9 -ye-y 9 + ye) - £4(2/7 ~ 2/7) + 
£7(2/4 - 2/4) + (£6 - £9X2/1 - 2 /i) + Ai(x 6 - £ 9 ) + 
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-£7(2/7 - 2/7) + AX®! - £4) , 


ix 4 = £4(2/3 - 2/9 - 2/3 + £9) - £7(2/1 - 2 /i) + 
£1(2/7 - 2/7) + (£9 - £3)( 2/4 - 2/4) + 
~K X 7 + A!, (£9 - £ 3 ), 

ix 6 = £7(2/7 - 2/7) - £7(2 17 - 2/7) + £1(2/1 - 2 /i) 
—£1(2/1 - 2/1) + Aj(£i - £1) , 
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